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Abstract It is believed that soon after the Planck time, Einstein's general 
relativity theory should be corrected to an effective quadratic theory. Nu- 
merical solutions for the anisotropic generalization of the Friedmann "open" 
model H 3 for this effective gravity are given. It must be emphasized that al- 
though numeric, these solutions are exact in the sense that they depend only 
on the precision of the machine. The solutions are identified asymptotically 
in a certain way. It is found solutions which asymptote de Sitter space, Rie- 
mann flat space and a singularity. The question of isotropisation of an initially 
anisotropic Universe is of great importance in the context of cosmology. Al- 
though isotropisation is not directly discussed in this present work, we show 
that sufficiently small anisotropies, do not increase indefinitely according to 
particular quadratic gravity theories. It can be understood as weak isotropisa- 
tion, and we stress that this result is strongly dependent on initial conditions. 
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1 Introduction 

The semi-classical theory considers the back reaction of quantum fields in a 
classical geometric background. It began about forty years ago with De Witt 
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PP, and since then, its consequences and applications are still under research, 
see for example [2]. 

Different from the usual Einstein-Hilbert action, the one loop effective grav- 
itational action surmounts to quadratic theories in curvature, see for exam- 
ple [T], (3]. It is the gravitational version of the Heisenberg-Euler electromag- 
netism. As it is well known, vacuum polarization introduces non linear cor- 
rections into Maxwell electrodynamics @], first obtained by Heisenberg- Eulcr 

This quadratic gravity was previously studied by Starobinsky [5] , and more 
recently, also by others [7] . The higher order terms are in fact responsible for a 
slowly decaying dark energy. In [7] only the homogeneous and isotropic space 
time is studied. It is of interest for example in the context of the final stages 
of evaporation of black holes, inflationary theories [5], in the approach to the 
singularity [pj and also in a more theoretical context [ID] . 

The effective gravity was apparently first investigated in Tomita's article 
[TT] for general Bianchi I spaces. They found that the presence of anisotropy 
contributes to the formation of the singularity. Gurovich and Starobinsky stud- 
ied the fourth-order f{R), in which the anisotropic part of the metric was 
integrated explicitly. Berkin's work shows that a quadratic Weyl theory is less 
stable than a quadratic Riemann scalar R 2 . In particular, in the very interest- 
ing article [12] , Barrow and Hervik addressed the anisotropic cases of Bianchi 
I and 77. The most interesting results of Barrow and Hervik are the exact 
solutions for quadratic theories of the same type investigated in this present 
article. Instead of the metric, the field equations in [12] are written in a differ- 
ent set of variables which by now is a well known procedure used in the context 
of dynamical systems in cosmology [T3]. We intend to switch to these set of 
variables in future works. Homogenous solutions in the context of quadratic 
gravity was also addressed by |14j . Schmidt did a review of higher order gravity 
theories in connection to cosmology [15]. 

Also in the context of quadratic theories we have the interesting alternative 
of Gauss-Bonnet type F(G) by Odintsov, Nojiri and collaborators [16], for 
recent reviews see [17] . When the action depends on a arbitrary function of 
the Gauss-Bonnet term, it is not a top invariant and a consistent dynamic 
follows from it. Theories of R 2 type are also investigated, for example in [18] , 

In our case the Gauss-Bonnet term is understood as a surface term and 
discarded. In fact, using the Gauss-Bonnet term and the relation between the 
Weyl tensor and Riemann tensor it is possible to demonstrate the equiva- 
lence between the particular theory under investigation in this work and other 
quadratic theories of gravity. 

We have previously numerically investigated Bianchi / type solutions [15] . 
Soon after this result we also investigated the stability of the this particular 
Bianchi / case [20J. It turns out that for zero cosmological constant A = 0, 
Minkowski geometry is structurally stable in the sense that there is a basin of 
attraction to Minkowski solution. And for A > de Sitter type geometry is 
structurally stable, also in the same sense. Thus, soon after the Planck era, the 
effective one loop quadratic gravity predicts for the very particular Bianchi I 
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models, that there is a basin of attraction to Minkowski space for A = and 
a basin of attraction to de Sitter space for A > 0, see [20] . 

Within Einstein's general relativity theory, Wainwright et al [2T] discovered 
that almost any initial condition for Bianchi VII a evolve to an exact plane 
wave solution. For a very good and recent review on homogenous models in 
the context of Einstein's theory, see Ellis [2"2"] . 

In this present work, the solutions for the Bianchi VII a spatial homogenous 
space for the effective quadratic gravity is addressed. This is the anisotropic 
generalization of the hyperbolic H 3 Friedmann "open" space. For any spatially 
homogenous space-time, the dynamical equations of motion, of any theory, re- 
sult in a non linear system of ordinary differential equations. For the particular 
quadratic gravity, the ordinary differential equation is of degree 4. 

The initial conditions chosen are very near to exact solutions of the quadratic 
theory, and also only classical vacuum source is considered, which we believe is 
the best description soon after the Planck era. The solutions are obtained nu- 
merically and they are understood asymptotically in a certain sense described 
in the text. It is shown that depending on the initial conditions and param- 
eters, the de Sitter solution, or Riemann flat solution, or a singularity type 
solution is obtained asymptotically. 

The paper is organized as follows. In section 2 a brief review of the coor- 
dinate free method is used to obtain the field equations for the Bianchi VII a 
case, and also the Newman-Penrose coefficients. In Section 3 3 different initial 
conditions are analyzed: one that converges to de Sitter space, one that con- 
verges to Riemann flat space, and one that converges to a singularity. In the 
last section we present our conclusions. We believe the index notation should 
be clear from the text. 



2 Theory and Development 

For completeness, we give a brief overview about the tetrad basis, coordinate 
free description of spatially homogenous geometries, for further details see for 
example [25] . 

The right-invariant basis which describes the group of isometries G3 are 
the Killing vectors that follow from the relation 



where CLu are the structure constants of the group which satisfies the Jacobi 
identity 



(1) 



C \ab C c )d - °- 



(2) 



Since right and left actions commute, 



[Va,£b] = 0, 



(3) 
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Table 1 BVIIa basis and structure constants. 
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where {rj a } is the left-invariant vector basis. The right-invariant 1-form basis 
{uj a }, its dual, follows from the first Cartan equation 

du a = -r a h Au b =r a bc u b Au c , (4) 

where the connection 1-form is defined by r b = r bc oo c . Since 27"^ = C a bc , 
then 

^ = 2^A/. (5) 

The Bianchi classification of all possible three-dimension algebras, based 
on the structure constants, can be obtained by Schiicking's method. Once the 
structure constants are known, the vector and 1-form basis are obtained in the 
following way: the first-order differential equation §5§ is solved to determine 
the basis {w a } and the duality relation uj a rib — S a b gives the basis {rj a }- Since 
£ €o (wV) = (A." 6 )»?c + u b (£(aVc) = -Ce a (* 6 c ) = and C ia rj b = = 

(from ([3])), then the Killing vectors {£ a } are obtained by solving the following 
first-order differential equation 

C^^i^i+^lddx^O. (6) 

The structure constants and the basis for the type Bianchi VII a {A is the 
group parameter) are shown Table [I] 

The 1-form basis performs the correct projection to obtain a spatially ho- 
mogeneous four-space 

ds 2 = g a b(t)u a <E> w h = -dt 2 + g a p(t)u) a <g> uA (7) 
In this new basis, it can be seen that, indeed, £ a is a Killing vector: 

C^g = C^{gabOJ a ®uJ b ) = {C ic g ab )u a ®uj b +g ab {C^ a )®uj b +g ab uj a ®{C ic uj h ) = 0, 

(8) 

remind that g a b(t) is a function of t only. The partial derivatives g%jk are 
replaced by directional derivatives g ab \ c = g a b,iVc' which are projections of the 
derivatives in the basis {r] a ,uj a }. The covariant derivative written in tetrad 
indices is 

V a T b = T b + r b ca T c , (9) 
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where the torsion- free connection follows from V a {gbc 

uj b ® uj c ) = and Q, 



r a bc — ^ (9ab\c + 9ac\b ~ 9bc\a + G a bc ~ 

The Riemann tensor is obtained from 
which according to ([9| and ( 10 ) results in 



a 



cab 



Cbacj- 



n bed — 1 bd\c ~ 1 bc\d ' 



r% r ld 



pa pi / 
1 fd 1 be 



^ cd 1 bf- 



(10) 



(11) 



(12) 



The Ricci tensor and Ricci scalar are obtained as usual: R a b = R c ac b an d 

R /'";,• 

The field equations for the semiclassical theory are obtained performing 
metric variations in the gravitational Lagrangian |19) . |20) 



C 



-A + R + a[ R ab R 



3^ 



PR 2 



(13) 



where C q is the quantum part of Lagrangian and a and j3 are constants. For the 
spatially homogenous space they are described by the tensor E — E a i,uj a ® 0J b , 



U! 



4 _ 



dt. 



where 



E ab = G ab + \gabA - ( /3 



a H 
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ab 



T„ 
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(14) 



Gab — Rab — ^9abR, 



Hob " 2^ ab ^ Ca ^ cd ~~ ^Rab 



H a b " ^9abR 2 — 2RR ao — 2g ao DR + 2R. a b, 



g a bOR + R-ab — 2R cd R c bda 



and the connection used is (10), the Riemann tensor is (12) and □ = V a V a is 
D'Alembertian operator. 

The second quantization program of quadratic gravity is provided by Stelle 
[21]. According to [21], the linearized filed equations result in eight degrees 
of freedom. Two correspond to the familiar massless spin 2 graviton. Five 
correspond to a massive spin 2 particle with mass mi = 1/ ' \fa. The last degree 
of freedom correspond to a massive scalar particle with mass mo = 1/ \f— 6/3. 
In this work we choose a > and j3 < throughout. The presence of tachyons 
occur if a < or /3 > and in our case, would indicate the linear instability 
of the solution. The existence of oscillations in our solutions is consistent with 
a normal scalar and a tensor massive particles. 

The counterterms i? 2 , R a bR ab , A and R in (13) are precisely the ones 
necessary in order to obtain a finite vacuum expectation value of the energy 
momentum tensor, see for example |25j . A theory without these counterterms is 
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inconsistent from the point of view of the renormalization of the quantum field 
in C q . The renormalized vacuum expectation value of the energy momentum 
tensor is set to zero, which emphasizes the effects of a theory that should have 
been considered from the start. We are disregarding any classical source in 
this work. 

The finite contributions to the vacuum expectation value of the energy 
momentum tensor are known only for very particular situations, split rank 
spaces, for which the heat kernel can be obtained exactly [26], then the point 
splitting method gives (0|T a b|0) exactly. These solutions are known as self 
consistent, G a b — (0|T a b|0) and there are very few cases, for example [2"7j . 

As in any other metric theory, the covariant divergence of E a b must be zero 

V a E ab = 0. 

From V a E a b = it follows that if E44 = and E^a = initially, then they will 
remain zero fo any time, and act as constraints on the initial conditions. Con- 
sequently, these constraints are checked to test the accuracy of the numerical 
results, while E a p = represent the real dynamical equations of the problem 
(see for example [2S], p. 165). 

The homogeneous and anisotropic metric supposed throughout this work 
in the numerical results is the following way, 



g a b{t) 



( a?(t) 

a 2 (t) a 4 {t) 

a 4 (t) a 3 2 (i) 
\ 






-v 



(15) 



ds 2 



-dt 2 



g a p{t)uj a <g)uA 



(16) 



where g a p is the spatial part, and the Bianchi VII a 1-form basis is given in 
Table [l] The replacement of the above line element in (14 1 taking into account 
the appropriate connection and covariant derivative defined in ( 10 ) and ^ 
results in a non-linear fourth-order ordinary differential equation system in 
the functions a,i(t), i — 1,2,3,4: 



dt 4 



di(t) = ft (o'j(t), a,j(t), a,j(t), aj(t)) 



(17) 



It is well known that the Bianchi VII a type contains the hyperbolic 
isotropic model. Replacing the Bianchi VII a basis Table [I] and instead of ( 1 



the isotropic metric g a p(t) = dia,g(a 2 (t),a 2 (t),a 2 (t)) into the line element 
( 16 1 results in 



ds 2 



-dt 2 + a 2 (t)[dx 2 + e 2Ax (dy 2 + dz 2 )} 



(18) 
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With the coordinate transformation 

ln(coshx 



sinh x cos ( 



y 



A 

sin 9 cos if 
A coth x — col 
sin 9 sin tp 



A coth x — cos 9 
the usual Friedmann line element is obtained, 



ds 2 = -dt 2 



A 2 



[d X 2 + sinh 2 X (d9 2 + sin 2 6 dip 2 )]. 



(19) 



First let us emphasize that every Einstein space, satisfying R a b = Ag a f,/2, 



is an exact solution of this one loop quadratic effective theory given in ( 14 1 
Note the particular case when the constant A = 0: vacuum solutions of Ein 



stein's equations are also exact solutions of (14 1. So there is the following 



exact solution of (14) for Bianchi VII a with the diagonal metric g a p(t) 
diag(af(t),a 2 2 (t),a7(i)) 



ai(i) = A 



a 2 (t)=AQ 



a 3 (t) = AG 




(20) 



where C\ and Ci are integration constants and A is the group parameter. This 
solution is in fact the de Sitter solution modulo a coordinate transformation. 
The particular case when A — > 0, [29] . 

ai(t) = At + d, 
a 2 (t) = C 2 (At + d), 
as (*) = C 2 (At + d), 
A^O, 

which is equivalent to Minkowski space since the Riemann tensor is identically 
zero, R% cd — 0. Morever, it is in fact Milne space 

ds 2 = -dt 2 + t 2 [d X 2 + sinh 2 X (d9 2 + sin 2 9dtp 2 )}, 

modulo a coordinate transformation. 

The intention is to characterize the obtained numerical solutions in some 



way. The Weyl tensor follows from the Riemann tensor ( 12 1 as 

Cabcd — Rabcd — x(.9ca-R&d + 9dbRca ~ 9cbRda ~ QdaRcb) + -^R{9ca9db ~ ffcfeffda) 
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A complex null basis can be denned 



k k a k t a k t a / l a / t a I t a t t a t t a 0, 

t a i a = -k a l a = 1, (22) 



with the corresponding null metric 
9ab = 9abA a B b 



I -1 0\ 
-10 00 
1 
\ 10/ 



(23) 



where A a and B are the null vectors in (22). The Newman-Penrosc complex 



coefficients are in fact the tetrad components of the Weyl tensor 

V^o = C abcd k a t b k c t d , 
= C abcd k a l b k c t d , 

^2 = C abcd k a t b t C l d , 
V>3 = C abcd k a l b i C l d , 
1p4 — C abcd t a l b t C l d . 

Also the Ricci tensor can be projected to the null tetrad, resulting into the 
independent components, 

Rkk = R a bk a k h 

Rkl = Rabk a l b 
Rll — Rabl a l b 
Rkt — Rabk a t b 
Rlt = Rabl a t b 
Rtt — Rabt a t b 

Rti = Rabt a i b , (24) 

the first 3 and the last one are real, and the remaining ones are arbitrary 
complex numbers comprising of course 10 independent components. 

The Newman-Penrose coefficients are related to the Petrov classification 
as shown for example in [23] and [28]. When all the i/j's are zero it's a Petrov 
type O, which characterizes a conformal Minkowski space and the Weyl tensor 
vanishes. An exactly conformally flat solution is either a generalized Friedmann 
solution or an interior Schwarzschild solution as it can be seen for example 
in [23] pg. 413. Minkowski and de Sitter geometries, are particular cases of 
conformally flat solutions. 

Before getting into the numeric result, we call the attention to the book 
[T5] . pgs. 62- 64. As it is well written there, there is not a complete statement 
as to what constitutes a minimal set for ensuring that a cosmological model is 
close to Friedmann-Lemaitre model. This book is based on an article 1301 for 
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which some assumptions are made: i) Einstein's equations are satisfied, ii) the 
source is a mixture of radiation and dust. 

We are not concerned if the two above conditions are satisfied in this 
present work. So we do not expect that their result, albeit being very in- 
teresting, should be verified in the particular context being discussed here. 

In the following we will present a particular example. Consider the Bianchi 
VIIa metric written in the appropriate base given in Table [l] 



( a(t) 2 \ 

b{t) 2 

b{t) 2 

\ -1/ 



(25) 



The time like vector u a — (0,0,0, 1) is geodesic and orthogonal to the group 
orbit. The magnetic part if a & — 0, and the electric E a f, part of the Weyl tensor 
is 

/2V> 2 0\ 
-ip 2 o 
-ip 2 
\ 0/ 



E b n = 



where ip2 is the above defined Newman-Penrosc coefficient 



^2 



1 / -b(t) 2 ii{t) + b(t)b(t)a(t) + b(t)a(t)b(t) - a{t)b{t) 2 
6 \ a(t)b(t) 2 



(26) 



Also regarding (25), the expansion H — 1/3(9, — V c u c is 



ld(i) 2bify 
3 a(t) 3 b{t) ' 



(27) 



Now consider 



o(t) 
b(t) = 



tA 

Z 1 -t-V3 sin(t) 
tA 

l-t-^coait)' 



(28) 



From ( 28 1 it is easily seen that ( 25 1 as t — >• oo is 



lim g a b 

t— ^oo 



(t 2 A 2 ^ 
t 2 A 2 
t 2 A 2 
\ -1/ 



(29) 



which according to (19 1 and (18) is Milne's space which is identical to the 
Minkowski solution. 
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Considering ( 28 ) and ( 25 1 , now we can explicitly obtain the asymptotic 



t — >• oo limits of the Newman-Penrose coefficient ( 26 ) and the expansion ( 27 ) 
to be 



t— > oo 



lim V2 = --£~ 1/4 cos(» 



f 



lim E ab E ab = -r L/2 cos(ty 

t—t-OO Q 

lim H = -^ 1/4 sin(i) 
t->oo 3 



(30) 



This example shows that Definition 2.f of [T3] 



^H ab H ab /H 2 « e 
VE ab E" b /H 2 « e 



(31) 



< e << f is not necessary for the asymptotic approach to a Friedmann- 
Lemaitre model. Of course the authors of [13] were aware of this fact as it can 
be read on pg. 64 where they carefully write that: if Definition 2.1 is satisfied, 
then the metric can be locally written in an almost-RW form. 



As a final word we will use a different criteria, namely, inspired in ( 30 ) , the 
asymptotic tetrad components of the Weyl and Ricci tensor 



lim Cabcd « e 

t— >oo 

lim R ab 

t— VOO 



(32) 



< e « 1 which is less restrictive than (31). We emphasize that the asymp- 



totic behavior of the tetrad components ( 32 ) is a weaker criteria than the one 
given in Definition 2.1 of [13] . 

For instance, the Kasner solution [3T] and the exact Bianchi VII a gravity 
wave [35] obey (32 1 with R ab = 0: on the other hand these two exact solutions 



do not satisfy (31 1. So in principle our criteria would not distinguish between 
Kasner exact anisotropic solution, and for example, the limit as t — > 00 given 



occurs in the strong sense. 



in ( 29 ) , with ( 28 ) and ( 25 1 , which is clearly a space for which isotropisation 



3 Numerical Solutions 

3.1 Asymptotically de Sitter Solution 



We choose an initial condition near the exact de Sitter ( 20 1 solution with t = 1 
and Ci = 0, C 2 = 20, A = 0.1, a = 1, = -5.0 and A = 0.02. The only non 
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null coefficients consistent with the E44 = and £41 = constraints are 



ai(t) = ^/| sinh f y ^ j - 2.8 x 10~ 3 



a 2 (t) = C 2 A J jsinhU 



03 (*) = C 2 W ^ sinh ( y 



o'i(t) = A cosh [ y — i 



a'2(t) = C 2 Acosh ( y -t 



a - 3(t) = C^Acosh 



di(t) = Ay ^ sinh 




d 2 (t) = C 2 A^smhU^t 

o 3 (t) = ^Ay^sinh fy^t 
a'i(i) = -1.14153106 x 1(T 2 
o 2 (t) = C 2 A^ cosh f\Zf*J 

a 3 (t) = -4.65026064 x HT 1 , (33) 

we stress that numerically the initial condition 04 = but it rapidly increases. 
The metric does not remains diagonal. 

The time evolution of the non zero Newman-Penrose coefficients, ipO'^i are 
shown in Fig. [l] For the line element chosen ipi = ^3 = 0. It can be seen that 
asymptotically the Newman-Penrose coefficients all vanish so the Weyl tensor 
is asymptotically zero. 

The time evolution of the non zero components of the Ricci tensor accord- 
ing to the null base (22) are shown in Fig. [2] For the particular metric chosen, 
(15) -(16), Rkt — Rit — 0. It is also shown in Fig. [2] that asymptotically, 
the Ricci tensor is a constant proportional to the null metric ( 23 ) . The only 
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Fig. 1 Numerical evolution of the real and imaginary parts of the Newman-Penrose coef- 
ficients ipQ - ip^. For the initial condition given in the text. The numerical integration was 
done until the proper time t = 100. 



consistent result with the held equation ( 14 ) is 
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We have numerically checked up to t — 500 in proper time that R a b — A/2g a b 
with one part in 10 8 . Also the constraints E44 = and E41 = where numeri- 
cally verified in Figure [3] which is a strong indication that the numerical result 
should be trusted. 
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Fig. 2 Numerical evolution of the non zero real and imaginary parts of the null tetrad ( \'2'2\ 
components of the Ricci tensor | |24[ |. The integration was done in the proper time. It can be 
seen that the Ricci tensor is asymptotically proportional to the null metric 1231. 
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Fig. 3 Numerical evolution of the constraints £44 and .E41 which should be zero identically, 
indicating that the numerical solution is accurate. The integration was done for the initial 
condition and parameters specified on the text and until the proper time 200. 
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Table 2 The only non null initial conditions together with j34[ l, consistent with the E44 
and E41 constraints 

di(0) d 2 (0) 03(0) 01 a 3 

A C 2 A C 2 A 2.45037249 x 1(T 2 9.82173813 x 10" 2 , 



3.2 Asymptotically Riemann Flat Solution 

In the numerical solutions the following values for the parameters C\ = 0.3, 
C* 2 = 2.0, A = 0.1, a = 1000.0, (i = -5.0, A = 0, where chosen, and initial 



conditions (t = 0) near Minkowski solution (21 1 



ai (t) = At + d = Ci 

a 2 (t) = C 2 (At + Ci) - 0.0015 = C 2 Ci - 0.0015 

a 3 (t) = C 2 {At + C l ) +4.0 x 10~ 4 = C 2 Ci + 4 x 10~ 4 

a A {t) = 0.0, (34) 

we stress that numerically the initial condition 04 = but it rapidly increases. 
The metric does not remains diagonal. 

The only non null initial conditions for the derivatives are given in following 
Table [2J and as before, the numerical values of d 3 ai/dt 3 and (Pa^/dt 3 are 
consistent with the two constraints S44 = and E41 = 0. 

The time evolution of the non zero Newman-Penrose coefficients, ^o~ip4 
are shown in Fig. [4] For the line element chosen = ips = 0. The time 
evolution of the non zero components of the Ricci tensor according to the null 
base (22) are shown in Fig. [5] For the particular metric chosen, (15) -(16), 



Rm = Ru — 0. 
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Fig. 4 Numerical evolution of the real and imaginary parts of the Newman-Penrose coef- 
ficients ipQ - ip^. For the initial condition given in the text. The numerical integration was 
done until the proper time t = 60. 



The time evolution of the constrains E44 and E41 are shown in Fig. [6| As 
the constraints are numerically verified, we believe the numerical solution is 
accurate. 

Since all the ip — > 0, Fig. [4j the Weyl tensor C£ cd — > 0, and since the Ricci 
tensor i? & — > 0, Fig.[5j then the Riemann tensor R a bcd — > 0. 

The numerical integrations shown in Fig. [4] and Fig. [3] were carried up to 
times t = 1 x 10 6 showing that null tetrad components of the Riemann tensor 
Rabcd — > in one part in 1 x 10 6 and smaller values. It is in this sense that 
the solution is understood to asymptotically approach Minkowski space. 
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Fig. 5 Numerical evolution of the non zero real and imaginary parts of the null tetrad i \'2'2\ 
components of the Ricci tensor \24\ . The numerical integration was done until the proper 
time t = 200. 
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Fig. 6 Numerical evolution of the constraints E44 and E41 which should be zero identically, 
indicating that the numerical solution is accurate. The integration was done for the initial 
condition and parameters specified on the text and until the proper time t = 200. 
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Table 3 The only non null initial conditions together with j35[ l, consistent with the E44 
and E41 constraints 

di(0) d 2 (0) 03(0) 'a-i «3 

A C 2 A C 2 A -5.03312185 x 10 2 -2.01311859 x 10'^, 



3.3 Singularity 

Using exactly the same parameters of the preceding section C\ = 0.3, C2 = 2.0, 
A = 0.1, a = 1000.0, (5 — —5.0, A — 0, a slightly different initial condition is 
chosen 

ai (t) = At + d = Ci 

oa(t) = C 2 {At + Ci) + 0.2 = C 2 Ci + 0.2 

a 3 (t) = C 2 {At + C l ) +4.0 x 10~ 4 = C 2 Ci + 4 x 10" 4 

o 4 (t) = 0.0, (35) 

we stress that numerically the initial condition 04 = but it rapidly increases. 
The metric does not remains diagonal. The only non null initial conditions for 
the derivatives are given in following Table [3] 

This initial condition evolves very fast to a singularity characterized by the 
increase of the curvature scalars R a bR ab and R a bcdR ahcd shown in Fig. [7J The 



0.02 [).[)4 0.0G 0.08 0.1 

t 

(a) (b) 

Fig. 7 a) The increase of the scalar curvature invariants R a t > R ab , R a bcdR abcd characterizing 
a singularity, b) The principle values of the metric showing a decrease of all the spatial 
directions, which we understand as a collapse. Where Ai = gn = af. A2 = [a| + &§ + v2]/2, 
A 3 = [a| + a\ - VA}/2 and A = a| - 2a§a§ + a| + 4a|. 




constraints which should be zero S44 and -E41 shown in Fig. [HJ prove that the 
numerical result is accurate. Note that as the singularity is approached, the 
numerical errors increase, which is expectable. Also, the numerical integration 
can be carried further and further, and the constraints are not really satisfied 
at all, and is not shown since we believe the result should not be trusted. 
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Fig. 8 The constraints .E44 and £41. Note that as the singularity is approached, the nu- 
merical errors increase, which is expectable. 



4 Conclusions 



In this present work, the numerical solutions of Bianchi VII a type are analyzed 
in the effective gravity context. This is the anisotropic generalization of the 
hyperbolic H 3 Friedmann "open" space. The quadratic gravity is the result of 
vacuum polarization counter terms which must be introduced into Einstein's 
theory of gravitation. It should be the most natural theory just after the Planck 
era. It must be emphasized that although numeric the solutions are exact in 
the sense that they depend only on the precision of the machine. 

Within the cosmological context it is of great interest the concept of isotropi- 
sation, in that an initially irregular Universe turns out less irregular by phys- 
ical processes. Or at least, the irregularities do not grow indefinitely. The 
supposition of a homogenous Universe is a very strong one, and allowing for 
cosmological anisotropy provides a less restricted initial condition than the 
isotropic Friedmann "open" H 3 model. As a less restricted initial condition 
makes it a more probable one for the creation of the Universe, whatever this 
means. We can conclude that weak isotropisation occurs in the sense that the 
anisotropics either vanish or at least do not grow indefinitely, depending on 
the initial condition, of course. According to quadratic gravity, not all initial 
conditions evolve to reasonable Universes. There are initial conditions that 
evolve to a singularity. 

Apparently, the question of stability of de Sitter solution in quadratic grav- 
ity was first obtained by Miiller, Schmidt and Starobinsky [55] . 

We have previously, numerically investigated Bianchi / type solutions [TI5] • 
We found that there is a basin of attraction to Minkowski space for A = 
and a basin of attraction to de Sitter space for A > 0, sec 20J. In this sense, 
Minkowski and de Sitter solutions are structurally stable according to the 
effective gravity, for the particular Bianchi / models we analyzed. Also in the 
Bianchi / case, there is a basin of attraction to a singularity. 

The numerical solutions are given in section 3. In section 3.1 the initial 
condition is chosen near the de Sitter exact solution. In section 3.2 and 3.3 the 
solution is chosen near Minkowski exact solution. Albeit initially the metric 
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is supposed diagonal in all cases studied, it gets non diagonal dynamically 



through the 04 (i) component, see eq. (151 



The numerical solutions are characterized asymptotically according to the 
null tetrad components of the Weyl tensor and Ricci tensor. We emphasize 
that this is not such a strong criteria as the one discussed in the end of section 
2. According to this criteria, for instance, Kasner's solution [3T] asymptotes 
Minkowski space. 

In this sense, for the first solution under consideration in section 3.1, the 
Weyl tensor is zero and the Ricci tensor is proportional to the metric R a b = 
A/2g a b asymptotically and we identify this solution with de Sitter space. In 
section 3.2 the solution asymptotes a Riemann flat space in sense that the Weyl 
tenosr and Ricci tensor vanish asymptotically. The convergence to Riemann 
flat space is not so fast and we shall not address this question in this present 
work. In section 3.3 a solution that converges to a singularity is presented. 

The existence of initial conditions that converge to a singularity, for in- 
stance the one given in section 3.3, is a strong indication that this quadratic 
gravity is certainly not a complete theory. 

In all the cases we considered, the numerical behavior described in the text 
was checked for much larger times that the ones plotted, and we believe that 
asymptotical interpretation is the correct one. 

Since the initial conditions are near de Sitter and Minkowski solutions, we 
can speculate that also for Bianchi VIIa solutions, de Sitter and Minkowski 



space should be structurally stable in the weak sense specified near eq. (32 1, 
according to the effective gravity, such that there should be basins of attraction 
to these solutions. We intend to address this stability issue in a future works. 
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